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NUCLEAR MAGNETIC RESONANCE IMAGING 
USING PHASE ENCODING 

BACKGROUND OF THE INVENTION 

5 

1. Field of the Invention 

This invention relates to the field of wellbore logging and, more particularly, to a method 
and apparatus for determining nuclear magnetic resonance logging characteristics of earth 
formations surrounding a wellbore, as a function of angular position about the borehole, either 
|g during the drilling of the wellbore or after drilling. 

tega±. 

iff 

jV 2. Description of Related Art 

ff . 

Hydrocarbon fluids, such as oil and natural gas, are obtained from a subterranean 

L geologic formation, referred to as a reservoir, by drilling a wellbore that penetrates the 

P, hydrocarbon-bearing formation. An understanding of the reservoir physical properties, often 

GI referred to as formation evaluation, is needed to determine the well's productive capacity, 

ru 

recoverable reserves, size and type of production equipment needed, and many other issues 
relating to the well's drilling, completion and production. Specific reservoir properties that are 
desired include, for example, porosity, permeability, and water saturations. % 

Electrical "logging" dates back to 1912, when Conrad Schlumberger began studying the 
20 problem of exploring the underground by means of surface electrical measurements. Since that 
time, various types of well logging techniques have been developed, such as acoustic, 
temperature, resistivity, nuclear, and gamma-ray measurement techniques. 
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Nuclear magnetic resonance (NMR) imaging involves the physical principle that various 
nuclei will precess at different frequencies in an imposed magnetic field. Many nuclei have a 
magnetic moment and behave much like a spinning bar magnet. This spinning motion is often 
referred to in the NMR terminology as the particle spin and is illustrated in Figure 27. Hydrogen 
5 nuclei, consisting of a single proton, has a relatively large magnetic moment and is found in both 
water and hydrocarbons that are located within a reservoir matrix. An externally applied 
magnetic field can interact with the spinning hydrogen protons and can produce measurable 
effects. An NMR logging tool can be designed to operate at the magnetic resonant frequency of 

j£ hydrogen, thereby allowing the tool to alter and detect the responses of the hydrogen protons 

M 

!|| within the region of investigation. By altering and detecting responses, the tool can obtain 

IssS: . 

i: information relating to the water and hydrocarbons within the reservoir. 

c | 

■SB- 

A static magnetic field is generated by the tool to initially align the hydrogen protons in 
y, the formation fluids. An oscillating radio-frequency magnetic field is generated by the tool to 
St alter the hydrogen protons alignment and tip the protons in a transverse plane. The tipped 
II protons move in a precessional motion around the initial alignment position in a manner similar 
to a spinning top that precesses in the Earth's magnetic field, as illustrated in Figure 27. Various 
NMR measurements of these and other related effects can provide an indication of the amount of 
total fluid contained within the formation, and can be used to indicate the identity of the fluid, 
whether water, gas or oil. The measurements can also provide indications on the pore and grain 
20 size distribution of the formation matrix and whether the fluids are bound within the formation 
matrix or are capable of movement, and therefore, potentially producible. 

One approach to obtaining nuclear magnetic resonance measurements involves inserting a 



Page 3 of 53 



PATENT 

NMR tool within the wellbore and applying a locally generated static magnetic field B 0 , which 
can be produced by one or more electromagnets or permanent magnets. The spins of the 
hydrogen protons within the formation matrix near the tool are aligned with the applied field Bo 5 
generating a net nuclear magnetization as the spinning hydrogen protons precess about the 
imposed magnetic field B 0 , as illustrated in Figure 28. Nuclear spins of the hydrogen protons 
align with the applied field B 0 , generating a net nuclear magnetization. Applying an RF field, B h 
perpendicular to the static field B 0 , as illustrated in Figure 29, can change the angle between the 
nuclear magnetization and the applied field B 0 . The frequency of the RF field should be equal to 
the Larmor frequency given by ©o = Y B o where y is the gyromagnetic ratio. After application of an 
RF pulse, the magnetization begins to precess around B 0 and produces a detectable signal in the 
antenna. As the protons precess about the static field B 0j they gradually lose synchronization 
with each other, as illustrated in Figure 30. This loss of synchronization causes the magnetic 
field in the transverse plane to decay. Phase encoding is caused by inhomogeneities in the static 
magnetic field and by molecular interactions. The signals can be analyzed to detect nuclear 
magnetic resonance properties of the formation and provide information relating to porosity, free 
fluid ratio, permeability, and other properties of the formation. See U.S. Pat. Nos. 4,717,878 
issued to Taicher et al. and 5,055,787 issued to Kleinberg et al. 

Another approach to obtaining nuclear magnetic resonance measurements employs a 
locally generated static magnetic field, B 0 , which may be produced by one or more permanent 
magnets or electromagnets, and an azimuthally-oriented oscillating magnetic field, Bi, which 
may be produced by one or more RF antenna segments that transmit and/or receive from different 
circumferential sectors of the logging device. See U.S. Pat. Nos. 5,977,768 and 6,255,817 
assigned to Schlumberger Technology Corporation. 

Page 4 of 53 



PATENT 

U. S. Pat. No. 5,796,252 issued to Kleinberg et al. describes a nuclear magnetic logging 
device that includes permanent magnets, an RF antenna, and a coil for generating a magnetic 
field gradient. The technique described in the '252 patent utilizes pulsed magnetic field gradients 
to obtain information regarding diffusion properties of the formation fluids. If internal gradients 
5 are present in the formation, a pulse sequence is applied to reduce or substantially eliminate the 
effect of internal gradients in the formation. The '252 patent does not identify a method for using 
the coil to obtain an azimuthal NMR measurement. 

U. S. Pat. No. 5,212,447 issued to Zvi Paltiel describes a nuclear magnetic logging device 

0 

f| that includes permanent magnets and an RF antenna coil. The '447 patent requires a magnetic 
field gradient coil to determine a diffusion coefficient, i.e., the rate at which molecules of a 

kvj ■ 

.£1 material randomly travel within the bulk of the same material. The ' 447 patent employs the 
CI diffusion coefficient to determine at least one of the following petrophysical parameters: 
water/hydrocarbon discrimination, water and hydrocarbon saturation levels, permeability, pore 

JUS*- 

fjj size and pore size distribution, oil viscosity, a measure of the average increase in electrical 
15 resistance due to the formation tortuosity, and q-space imaging of the formation. The '447 patent 
does not identify a method for using the coil to obtain an azimuthal NMR measurement. 

U.S. Pat. No. 6,326,784 assigned to Schlumberger Technology Corporation, discloses a 
means to obtain azimuthal NMR measurements using one or more gradient coils and an axi- 
symmetric antenna. In this approach, a pulse sequence and a firing of a gradient coil is used in 
20 conjunction to obtain an azimuthal image. The resolution of the azimuthal image depends on the 
number of gradient coils used and the angular coverage of each gradient coil. The gradient coils 
are positioned circumferentially and are separated by an angular distance, for example, three 
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gradient coils located around an NMR tool spaced 120° from each other. Each gradient coil is 
used to spoil or rotate the hydrogen proton spins within the formation matrix adjacent to the 
gradient coil, with negligent effects everywhere else. The NMR data obtained after the pulse 
sequence and firing of a gradient coil is used to obtain formation evaluation information of the 
reservoir adjacent to the gradient coil This process is repeated for the other gradient coils to 
obtain the azimuthal image. 

One method for acquiring azimuthal data is the binning of the measured data. In this 
acquisition scheme, a plurality of azimuthal bins is defined and each NMR measurement is added 
to the content of a buffer associated with the bin in which the measurement was taken. To obtain 
a good statistical quantity of data for each bin and to ensure that in each bin there are enough 
measurements, this scheme requires numerous individual measurements. For example, for one 
particular embodiment having a 90° resolution, seven bins are preferable, with ten measurements 
per bin, or 70 measurements per scan. Better resolution, (i.e., smaller angular distances between 
gradient coils), would need substantially more measurements and result in substantially more 
time required to obtain the quantity of data needed. 

There is a need for improved means to obtain NMR azimuthal images of a formation. 
There is also a need for greater azimuthal resolution of the formation properties, along with 
means to obtain this information with less acquisition time. 
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SUMMARY OF THE INVENTION 

One embodiment of the present invention is a method for nuclear magnetic resonance 
imaging of an investigation region of formation surrounding a wellbore. The method comprises 
the steps of applying a series of magnetic field gradients to phase encode nuclei spins within the 
5 investigation region, wherein the strength of the magnetic field gradient applied for each spin- 
echo signal is different from at least one previously applied magnetic field gradients within the 
series. Nuclear magnetic resonance signals are detected from the investigation region resulting 
from the series of magnetic field gradients. 

O Another embodiment of the invention is an apparatus for determining a nuclear magnetic 

^ resonance property in an investigation region of earth formations surrounding a wellbore. The 
,1 apparatus comprises a logging device moveable through the wellbore and a means in the logging 
B device for applying a static magnetic field circumferentially around the wellbore and into the 
M* investigation region. The apparatus further comprises an antenna means in the logging device for 
Cf applying an RF magnetic field circumferentially around the borehole and into the investigation 
15 region, whereby the antenna means induces a plurality of pulse echoes and spin-echo signals 
from selected nuclei of the formation. A gradient means in the logging device is capable of 
producing a gradient magnetic field circumferentially around the borehole and into the 
investigation region and is capable of producing different strength gradient magnetic fields for 
each of the plurality of pulse echoes, wherein the orientation of the gradient magnetic field and 
20 the static magnetic field effects on the selected nuclei vary depending on the azimuthal position 
around the wellbore in relation to the gradient means and means for detecting nuclear magnetic 
resonance signals from the investigation region. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

The invention may be understood by reference to the following description taken in 
conjunction with the accompanying drawings, in which like reference numerals identify like 
elements, and in which: 

Figure 1 shows a well logging apparatus in which embodiments of the invention can be 
practiced; 

Figure 2 illustrates a logging while drilling (LWD) apparatus in which embodiments of 
the invention can be practiced; 

Figure 3 illustrates a MWD tool that includes magnetometers H x and H y oriented along x 
H and y axes of the tool that can include embodiments of the present invention; 

° Figure 4 is a flow chart that describes the Quadrant/Coil Position Determination Program; 

O 

f* Figure 5 illustrates a nuclear magnetic resonance (NMR) logging-while-drilling tool in 

2 accordance with one embodiment of the invention; 

m 

Figure 6 is a diagram of the firing circuitry for use in conjunction with an embodiment of 
15 the invention having two gradient coils; 

Figure 7 illustrates a nuclear magnetic resonance (NMR) logging- while-drilling tool in 
accordance with one embodiment of the invention; 

Figure 8 illustrates a cross sectional view of a nuclear magnetic resonance (NMR) tool 
having the induced static magnetic field adjacent to the gradient coil, in accordance with one 
20 embodiment of the invention; 
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Figure 9 illustrates a cross sectional view of an embodiment of a nuclear magnetic 
resonance (NMR) tool and an induced magnetic field imposed to the gradient coil, in accordance 
with one embodiment of the invention; 

Figure 10 illustrates a cross sectional view of an embodiment of a nuclear magnetic 
5 resonance (NMR) tool and an induced magnetic field imposed to the gradient coil shown as an 
approximate average for each segment of the investigation region, in accordance with one 
embodiment of the invention; 

u Figure 1 1 illustrates a cross sectional view of an embodiment of a nuclear magnetic 

B resonance (NMR) tool and an induced magnetic field imposed to the gradient coil shown as an 

If I 

ft approximate average for each segment of the investigation region along with the induced static 
magnetic field, in accordance with one embodiment of the invention; 

Figure 12 illustrates the pulse sequence used to obtain the azimuthal resolution imaging 
for one embodiment of the present invention; 

HI Figure 13 is a table showing the NMR signals from various zones arranged in five bins, in 

1 5 accordance with one embodiment of the invention; 

Figure 14 illustrates a cross sectional view of an embodiment of a nuclear magnetic 
resonance (NMR) tool, in accordance with one embodiment of the invention; 

Figure 15 illustrates the pulse sequence used to obtain the azimuthal resolution imaging 

4 

for an alternative embodiment of the present invention; 

20 Figure 1 6 is a table showing NMR signal intensities resulting from the second gradient 

pulse sequence grouped in five bins, in accordance with one embodiment of the invention; 



yy 
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Figure 17 illustrates an embodiment of the present invention; 

Figure 18 illustrates the radial and angular directions of an embodiment of the present 
invention; 

Figure 19 graphically illustrates the amplitude of the radial component as a function of 
the azimuthal angle; 

Figure 20 illustrates an embodiment of the present invention; 

Figure 21 illustrates an alternate embodiment of the present invention; 

Figure 22 graphically illustrates the magnetic field profile versus radial distance of an 
embodiment of the present invention; 

Figure 23 illustrates the partitioning of a formation into radial segments in one 
embodiment of the present invention; 

Figure 24 illustrates an alternate embodiment of the present invention; 

Figure 25 graphically illustrates the magnetic field profile versus axial distance of an 
embodiment of the present invention; 

Figure 26 illustrates the partitioning of a formation into axial segments in one 
embodiment of the present invention; 

Figure 27 illustrates the concepts of nuclei spinning motion and precessional motion 
within a three-dimensional space; 

Figure 28 illustrates a directionally imposed static magnetic field within a three- 
dimensional space; 
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Figure 29 illustrates the concept of spin tipping and an applied transverse magnetic field; 

and 

Figure 30 illustrates the concept of dephasing of nuclei precessing around a static 
magnetic field. 

5 

DESCRIPTION OF ILLUSTRATIVE EMBODIMENTS 

Illustrative embodiments of the invention are described below. In the interest of clarity, 
q not all features of an actual implementation are described in this specification. It will of course 

III be appreciated that in the development of any such actual embodiment, numerous 

U.- 

# implementation-specific decisions must be made to achieve the developers' specific goals, such 

SI 

w as compliance with system-related and business-related constraints, which will vary from one 
2 implementation to another. Moreover, it will be appreciated that such a development effort 
CP might be complex and time-consuming, but would nevertheless be a routine undertaking for 
RJ those of ordinary skill in the art having the benefit of this disclosure. 

15 Figure 1 shows a well logging apparatus 10 in which embodiments of the invention can 

be practiced. The well logging apparatus 10 is disposed in a wellbore 12 drilled from the surface 
22 and it is capable of making measurements of properties of the various earth formations 13, 14, 
15, 16. The wellbore 12 is typically filled with a liquid 18 such as drilling mud or completion 
fluids. A region within a particular formation 15 exposed to nuclear magnetic resonance 

20 ("NMR") imaging is shown generally as 20 and it has a generally cylindrical shape extending 
within the formation and proceeding around the wellbore. 
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The well logging apparatus 10, which can include an NMR tool 30, can be lowered into 
the wellbore 12 by means of a wireline 24 that can provide electrical power and transmit signals 
to the surface 22. Surface equipment (not shown) can include a telemetry system for 
communicating with the tool string, recording instruments, computers, and the like. The well 
logging apparatus 10 is often centered within the wellbore 12 by means of a top centralizer 28 
and a bottom centralizer 29 attached to the well logging apparatus 10 at axially spaced locations. 
The centralizers 28, 29 can be of various types known in the art, such as bowsprings. 

Other well logging sensors may form part of the well logging apparatus 10. Additional 
logging sensors 32, 34 may be located above the NMR probe 30 while other logging sensors, 36, 
38 may be located below the NMR probe 30. The other sensors 32, 34, 36, 38, can be of types 
familiar to those skilled in the art and can include, but are not limited to, pressure and 
temperature sensors, gamma ray detectors, formation bulk density sensors or neutron porosity 
detectors. The locations of the other sensors 32, 34, 36, 38 are a matter of convenience and are 



2 not to be construed as a limitation on the invention. 



hi 

15 Referring to Figure 2, there is illustrated a logging while drilling (LWD) apparatus 40 in 

which embodiments of the invention can be practiced. A drill string 42 is disposed within the 
wellbore 12 and includes a drill bit 44 at its lower end. The drill string 42, and the drill bit 44 
can be rotated by a rotating table (not shown) that engages a kelly (not shown) at the upper end of 
the drill string 42. Alternatively, the drill string 42 may be rotated from the surface by a "top 

20 drive" type of drilling rig. Alternatively, the drill string 42 may comprise a downhole motor (not 
shown) that rotates the drill bit 44, without requiring the rotation of the drill string 42. Drilling 
fluid, often referred to as mud, is pumped downward through a channel in the center of drill 
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string 42. The drilling fluid exits the drill string 42 via ports in the drill bit 44 and then circulates 
upward in the region between the outside of the drill string 42 and the periphery of the wellbore 
12. The drilling fluid thereby carries formation cuttings to the surface of the earth. 

Tools designed for logging while drilling 50 (LWD), measurement while drilling 52 
5 (MWD), or a combination of both (LWD/MWD) can be connected to the drill string 42. A 
typical MWD tool 52 measures and/or computes the direction, inclination, and rotational 
orientation of the bottom hole assembly. An MWD tool useful with the subject invention is 
described, for example, in U. S. 5,473,158. The driving electronics module 54 and acquisition 
p and processor electronics module 56 are coupled to and obtain measurement information from 
jg the LWD tool 50. 

M Figure 3 illustrates a MWD tool 52 that includes magnetometers H x and H y (60, 62) 

oriented along x and y axes of the tool that can include or be used in conjunction with 
embodiments of the present invention. The MWD tool 52 can be used to gather alignment 
information relating the orientation of the present invention with surface orientations, such as the 
15 directional orientation of formation properties (e.g., a north-east by south-west aligned fractured 
zone through the formation). Such x and y axes are in the plane of a radial cross section of the 
tool. A z-axisofthe tool is oriented along its longitudinal axis. In a similar way, accelerometers 
G x and G y of the accelerometer package 64 (which also includes an accelerometer along the z- 
axis of the tool) are oriented along the x and y axes of the tool. A microcomputer 58 responds to 
20 H x and H y signals and G x and G y signals to constantly determine an angle <|> between an H' 
vector and the G f vector, in the cross sectional plane of the MWD tool 52. The H' vector 
represents that portion of a vector pointed to earth's magnetic north pole that is projected onto 



Pt 
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the x-y plane of MWD tool 52. The G' vector represents the down component in the cross 
sectional plane of tool 52 of the earth's gravity vector. A signal representative of such angle <j> is 
constantly communicated to the downhole computer 58 (which includes a Quadrant/Coil Position 
Determination program). The MWD tool 52 can be used in conjunction with the present 
5 invention to distinguish the azimuthal location and orientation of angular segments of the 
formation around the wellbore. Various embodiments of this type of MWD tool can be utilized 
in rotational drilling operations, in non-rotational drilling operations in which the drill string 
moves primarily longitudinally through the wellbore, and also in wireline conveyed logging 
B applications. 

it 

iffl Figure 4 is a flow chart that describes the Quadrant/Coil Position Determination Program 

'# 100. As explained above, an angle <|> is constantly computed between the H' vector (a constantly 

□ 

y[ directed vector in the x-y plane for a vector directed to earth's magnetic pole) and a G' vector ( a 

jjl constantly directed down vector in the x-y plane of a vector directed to the earth's gravitational 

LI 

W center). As the tool rotates in the wellbore, the x and y axes of the device rotate at the angular 
15 speed of the drill string so the x and y components of the H' vector and G' vector are constantly 
changing with time. Further, as the device rotates in the wellbore, an angle 0(t) is constantly 
formed between the tool x-axis and such H' vector. The angle 8(t) is determined from the H x 
and H y signals from magnetometers 60 and 62 and the angle varies with time because it is 
measured from the x-axis of the MWD tool 42 (and the LWD tool 50) to the H' vector. 

20 At step 102, the down vector angle, ZD(t), is determined in Quadrant/Coil Position 

Determination program 100, according to the following relationship, as a function of the x and y 
axes and time: 
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G(t) =cos" 1 



lH x {t) 2 + H y {tf) 



The angle of the down vector is determined in the program as ZD(t) = 0(t) - <|>. 

At step 104, four quadrants may be defined by angular ranges about the periphery of the 

tool: 

5 Q B o T (t) = ZD(t) -45° to ZD(t) +45° 

Q LE Fr(t) = ZD(t) +45° to ZD(t) +135° 

Q 

8 QtopCO = ZD(t) +135° to ZD(t) +225° 

ss?r? 

5 QRiGHx(t) = ZD(t) +225° to ZD(t) -45°. 



The term "quadrant" is used to illustrate the invention where four 90° angular distance segments 
W) are defined around the 360° circumference of the MWD device or the LWD tool Other angular 
|f{ distance segments, sometimes referred to as zones or segments, may be defined having 
resolutions other than 90°. The smaller the angle of investigation, the greater the formation 
resolution is obtained. 

Figure 5 illustrates a nuclear magnetic resonance (NMR) logging-while-drilling tool 50 in 
15 accordance with one embodiment of the invention. The tool 50 has an axis 60 and can comprise 
a portion of a drill collar 62 5 which is substantially aligned with the axis of the wellbore. A static 
magnetic field is produced by tubular, axially polarized, permanent magnets 64, 66 that are 
mounted inside the drill collar 62. A channel 68 located inside the tool permits drilling mud to 
flow toward the drill bit. In the region between the magnets 64, 66, there is a recessed area 70. 
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An RF antenna 72 is provided in the recessed area 70, which can be used for detecting NMR 
signals. However, a separate antenna or receiver may be used to detect the signals. A non- 
conductive material 74 can be provided in the recessed area 70 beneath the antenna 72. The 
material 74 may be a ferrite to increase the efficiency of the antenna 72. Alternatively, the 
5 material 74 may comprise a plastic, rubber, or a reinforced epoxy composite material. 

Still referring to Figure 5, in order to obtain azimuthal NMR measurements, at least one 
gradient coil 76 can be arranged in the recessed area 70. The geometry of the gradient coil 76 
1 1 (e.g., a saddle coil) can produce a magnetic field around the tool 50, as illustrated in Figure 9. 
§ The gradient magnetic field will add to or subtract from the static magnetic field, depending on 



its azimuthal orientation to the gradient coil. It should be noted that the method and gradient 
coils of the subject invention can also be used with any tool that generates a rotationally 



L symmetric static magnetic field, for example, the tools disclosed in U. S. Pat. Nos. 5,757,186 
issued to Taicher et al., and 5,280,243 issued to Melvin Miller. 



H Figure 6 is a diagram of an embodiment of the circuitry included in the driving 

15 electronics module 54 for use in conjunction with an embodiment having two gradient coils 76 
and 94. An example of a tool having two gradient coils is illustrated in Figure 14. The driving 
electronics include a high voltage power supply 80 and capacitor 82. Switches Si and S 2 are 
under the control of timing control/coil selection hardware 84. The gradient coils 76 and 94 are 
coupled, via the switches, to the timing control/coil selection hardware 84. The pulse sequence 
20 created by the gradient coils can also comprise applying a fixed or variable wait time between the 
gradient pulse and the RF pulse. 
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Figure 7 illustrates a nuclear magnetic resonance (NMR) logging-while-drilling tool 50 in 
accordance with one embodiment of the invention. A static magnetic field 88 is produced by 
tubular, axially polarized, permanent magnets 64, 66 that are mounted inside the drill collar 62. 
The static magnetic field 88 generated by the magnets 64, 66 extends out radially from the tool 
50 into the formation that is axially adjacent to the gradient coil 76. The static magnetic field 88 
can also be developed using an electrical current flowing through coils within the tool 50. 

Figure 8 illustrates a cross sectional view of a nuclear magnetic resonance (NMR) tool 90 
having the induced static magnetic field 88. The induced static magnetic field 88 extends out 
radially into the formation that is axially adjacent to the gradient coil 76. 

Figure 9 illustrates a cross sectional view of an embodiment of a nuclear magnetic 
resonance (NMR) tool 90 wherein the geometry of the gradient coil (e.g., saddle coil) 76 is such 
that it produces a polarized magnetic field 92 around the tool 90. The magnetic field 92 extends 
from one end of the gradient coil 76a, proceeds circumferentially through the adjacent formations 
on either side of the tool 90, and returns to the other end of the gradient coil 76b. The gradient 
magnetic field 92 is oriented circumferentially into the investigation region relative to the 
wellbore. 

Figure 10 illustrates a cross sectional view of an embodiment of a nuclear magnetic 
resonance (NMR) tool 90 wherein the geometry of the gradient coil 76 produces a magnetic field 
92 around the tool. The formation around the tool is shown divided into eight zones, each zone 
representing an azimuthal angle of 45°. The dashed arrows show an approximation of the 
average orientation of the gradient coil magnetic field 92 within each zone. The magnetic field 
92 extends from one end of the gradient coil 76a, proceeds through the adjacent formations on 
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either side of the tool 90, and returns to the other end of the gradient coil 76b. The gradient 
magnetic field 92 exerts differing directional magnetic forces, depending on its azimuthal 
position relative to the gradient coil 76 for each zone. The gradient magnetic field amplitude is 
not linear, but varies relative to the cosine of the azimuthal angle from the gradient coil, therefore 
5 the dashed arrows should not be interpreted to indicate the actual direction and strength or the 
magnetic field 92, but only a generalized indication. 

Figure 1 1 illustrates a cross sectional view of an embodiment of a nuclear magnetic 
resonance (NMR) tool 90 wherein the geometry of the gradient coil 76 produces a magnetic field 
Q 92 around the tool. The formation around the tool is shown divided into eight zones, each zone 

5 

j|j representing an azimuthal angle of 45°. The dashed arrows show an approximate average of the 
Si magnetic field orientation of the gradient coil magnetic field 92 on each zone. The magnetic 
s field 92 extends from one end of the gradient coil 76a, proceeds through the adjacent formations 
on either side of the tool 90, and returns to the other end of the gradient coil 76b. The gradient 
S[ magnetic field amplitude is not linear, but varies relative to the cosine of the azimuthal angle 
5 f5 from the gradient coil, therefore the dashed arrows should not be interpreted to indicate the actual 
direction and strength or the magnetic field 92, but only a generalized indication. The solid 
arrows show the static magnetic field 88. The gradient coil field 92 acts to increase the static 
magnetic field 88 by a maximum amount in zone 1, decrease the static magnetic field 88 by a 
maximum amount in zone 5, and have effects between these two extremes in the other zones. 
20 The eight zones are typically fixed in space relative to the gradient coil 76, but can be defined as 
a function of time if the tool is rotating, as in some embodiments of a logging while drilling 
(LWD) application, (see prior discussion involving Figures 3 and 4). 
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In the typical NMR measurement application, a pulse sequence is applied to the 
formation under investigation. In one embodiment of the invention a pulse sequence, such as the 
Carr-Purcell-Meiboom-Gill (CPMG) sequence, first applies an excitation pulse, a 90° pulse, to 
the formation that rotates the spins into the transverse plane (relative to the static magnetic field). 
5 After the spins are rotated by 90° and start to dephase, the carrier of the refocusing pulses, the 
180° pulses, is phase shifted relative to the carrier of the 90° pulse sequence according to the 
following relationship: 



y t 

H 5 where the bracketed expression is repeated for n=l, 2,...N, where N is the number of echoes 

£: 

M collected in a single CPMG sequence and the echo spacing is t echB = 2t cp = t m +t x +t 2 . The term 

P 90 °±x denotes an RF pulse that causes the spins to rotate by a 90° angle about the +x-axis, as 
E commonly defined in the rotating frame of magnetic resonance measurements (phase alternated). 
5 The time between application of the 90° pulse and the 180° pulse, t 0 , is r less than t cp , half the 

echo spacing. The CPMG sequence enables acquisition of a symmetric measurement (i.e., a 
15 measurement without using the gradient coils). The exact timing parameters, t 0 , t u and t 2 , 

depend on various factors (e.g., the shape of the applied pulses). 



Figure 12 illustrates the pulse sequence used to obtain the azimuthal resolution imaging 
for one embodiment of the present invention. Referring to Figures 8 - 1 1, a current pulse applied 
to the gradient coil 76 generates an additional magnetic field 92, altering the static magnetic field 
20 88 within the formation matrix, depending on its position relative to the gradient coil 76. The 
gradient coil current amplitude is different for each consecutive echo in the CPMG pulse 
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sequence, so that at each consecutive echo, the nuclear spin alignments within the formation 
matrix are rotated by a different angle. Within the present application the phrases of "different 
gradient current amplitude" or "altering the gradient amplitude" and the like can include 
incrementing, decrementing or other means of altering the amplitude such that it is not the same 
5 amplitude as at least one previous gradient coil pulses, thus imparting distinct pulse amplitudes 
on some or all of the altered echoes. The currents shown in the figures illustrate increasing 
current amplitudes, but other patterns of alteration can also be used. The amplitude of the 
gradient pulse that provides the greatest phase change has to be at least strong enough to impose 

|s«Ss 

p £ ... 

X an altered phase difference between adjacent formation segments (areas of investigation or 

li resolution) of at least 180 degrees, thereby completely dephasing neighboring formation 

£ 

K 4 segments. 

& 

p, Typically the gradient coil current amplitude is altered for each consecutive echo, but 

IsSSft 

I* keeping consistent gradient coil current amplitudes and increasing the incremental times between 
59 the pulses can obtain a similar effect. Generally shorter, more consistent time periods are 
15 preferred within the pulse sequence; therefore the current amplitude is generally the variable that 
is altered. The gradient coil current pulse is applied just before each echo. The gradient coil 
current is then reversed at the same amplitude after each echo so that phase encoding or 
rotational effects due to the gradient pulse are cancelled before the next 180° phase reversing 
pulse. 

20 Since the strength of the static magnetic field 88 is much greater than the gradient coil 

magnetic field 92, the only component of the gradient coil magnetic field 92 that has an effect on 
the nucleus spins is those components that are parallel or anti-parallel to the static magnetic field 
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88. Only the component of the gradient field 92 along the static magnetic field 88 orientation is 
responsible for the changes in the precession frequency (or phase) of the spins, i.e. that portion of 
the gradient coil magnetic field 92 that is parallel or anti-parallel in respect to the static magnetic 
field 88. 

5 For the embodiment of the invention using the pulse sequence as shown in Figure 12, the 

gradient coil pulse is applied to the five echoes in the CPMG sequence. The echo amplitude with 
respect to the gradient pulse amplitude is analyzed using Fourier transforms, which results in data 
*? that can be grouped within five bins representing NMR signal intensities. The NMR signal 
J intensities from the various bins correspond to either individual zones or combinations of the 

Ip zones of the formation located around the NMR tool. U.S. Pat. No. 6,326,784 assigned to 

St 

# Schlumberger Technology Corporation discloses one means of utilizing Fourier transforms of 

5 

O NMR signals that can be grouped within bins, the binned data representing groupings of NMR 

i t 

E signals obtained from one or more azimuthal zones around the NMR tool. 

rJ Figure 13 is a table showing how the NMR signal intensities from the five bins 

15 corresponding to zone 1, zones 2 and 8, zones 3 and 7, zones 4 and 6, and zone 5, respectively. 
Zones 1 and 5 can be individually distinguished since they each are in a bin by themselves, 
however the other six zones are in combinations within the bins and therefore cannot be directly 
distinguished. Since the gradient magnetic field amplitude is not linear, but varies relative to the 
cosine of the azimuthal angle from the gradient coil, the zones 1 through 8 will comprise 
20 differing azimuthal angles and will therefore provide data for differing angular segments of the 
formation. At this stage of the formation evaluation, the opposing zones 1 and 5, alternately 
referred to as top and bottom segments, can be quantified and are obtained using a single gradient 
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coil The time to obtain this data is the time required to generate and detect five CPMG echoes. 
To quantify other formation segments, in one embodiment of the invention, the NMR tool can be 
rotationally adjusted to a new azimuthal position, and the method described above repeated. 

Figure 14 illustrates a cross sectional view of an embodiment of a nuclear magnetic 
5 resonance (NMR) tool 90 where a second gradient coil 94 is oriented at 45 degrees (e.g. facing 
hypothetical zones 2 and 6) from the original gradient coil 76. The dashed arrows illustrate an 
approximate orientation of the magnetic field generated by the second gradient coil 94 on each 
zone. The solid arrows represent the static magnetic field 88. The various segments of the 
formation matrix are typically fixed in space relative to the gradient coil 94, but can be defined as 
f| a function of time if the tool is rotating, as in some embodiments of a LWD application. The 
*p CPMG sequence as shown in Figure 12 can be repeated for the second gradient coil 94, or the 
® second gradient coil firing can be performed within the same pulse sequence after the first 
rf gradient coil firing. 

CO Figure 15 illustrates the pulse sequence used to obtain the azimuthal resolution imaging 

¥5 for an alternative embodiment of the present invention. In this illustrative embodiment, the first 
gradient coil is fired with incremented amplitudes before and after echoes 1 through 5, while the 
second gradient coil is fired with incremented amplitudes before and after echoes 6 through 10. 
The gradient coil current amplitude is altered for each consecutive echo in the CPMG pulse 
sequence, so that at each consecutive echo, each segment is rotated by a different angle. 

20 The echo amplitudes with respect to the second series of gradient pulse amplitudes are 

analyzed using Fourier transforms, which results in an image with five bins. The NMR signal 
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intensities from the various bins correspond to either individual zones or combinations of the 
zones of the formation located around the NMR tool. 

Figure 16 is a table showing how the NMR signal intensities resulting from the second 
gradient pulse sequence can be grouped in five bins corresponding to zone 2, zones 1 and 3, 
zones 8 and 4, zones 7 and 5, and zone 6, respectively. Zones 2 and 6 can be individually 
distinguished since they each are in a bin by themselves. Again, it is important to note that the 
various zones represent differing azimuthal segments within the formation with differing 
azimuthal angles, for example, the area of zone 2 in Figure 16 will have a differing azimuthal 
angle of coverage than the area of zone 2 in Figure 13. It is possible for the formation coverage 
of zone 1 in Figure 13 to overlap the formation coverage of zone 2 in Figure 16, however, the 
ability to determine the NMR signal intensities of two adjacent zones, even with an overlap of 
areas, can be beneficial in providing formation data with greater azimuthal resolution. For 
example, the zone 1 as described above relating to Figures 12 and 13 provides data for an area 
with an average azimuthal angle in relation to a first gradient coil, while the zone 2 as described 
above relating to Figures 15 and 16 provides data for an area (part of which overlaps with zone 1 
from Figure 13) having a different average azimuthal angle in relation to the first gradient coil. 
The difference in the average azimuthal angles can be considered the formation resolution 
between these two sets of data, but does not provide this level of formation resolution around the 
entire wellbore circumference. 

Figure 17 illustrates an embodiment of the present invention in which azimuthal imaging 
can be obtained with two gradient coils 76, 94 that are located perpendicular to each other (i.e., 
90° spacing from each other). The analysis of the NMR signals relating to the first gradient coil 
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76 with five phase encoding gradients provides data representing a top zone 110 and a bottom 
zone 112, each having an azimuthal resolution covering approximately 106 . The analysis of the 
NMR signals relating to the second gradient coil 94 with five phase encoding gradients provides 
data representing a left zone 114 and a right zone 116, each having an azimuthal resolution 
covering approximately 106°. Other azimuthal resolutions can be obtained by altering the 
number of phase encoding gradients that are used, for example, the azimuthal resolution can be 
reduced from approximately 106° to approximately 83° by using eight phase encoding gradients 
rather than the five gradients discussed above. 

Azimuthal images like those that can be obtained from a second gradient coil located at 
45° resolution from the first gradient coil can also be obtained with two gradient coils oriented 
perpendicular (i.e., 90°) to each other as illustrated in Figure 17. The first phase encoding 
gradient coil 76 can be fired by itself, as described above, to produce a first set of binned results 
corresponding to the groupings of a first set of azimuthal zones. During the firing of the second 
gradient coil 94, the first gradient coil 76 can also be fired, both gradient coils fired at 
approximately the same time and having approximately the same amplitude. The result of the 
simultaneous firing of both gradient coils 76, 94 produces a gradient field that is comparable to a 
field that would have been created by a single gradient coil located between the two gradient 
coils (i.e., at an azimuthal position 45° from the first coil). The simultaneous firing of both 
gradient coils replicates the effects of a hypothetical single gradient coil located at 45° spacing 
from the first gradient coil and results in a second set of five different bins of data corresponding 
to the groupings of a second set of azimuthal zones. The orientation of the magnetic field 
generated by the simultaneous firing of both gradient coils can be varied by changing the relative 
current amplitudes of the two gradient coils as they are fired together, thereby obtaining 
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azimuthal resolutions approximately equivalent to a hypothetical gradient coil, located at an 
equivalent gradient coil phase angle, without physically altering the location of the gradient coils. 
The cumulative magnetic field gradient pulse effect of the two simultaneously fired gradient 
coils defines the equivalent gradient coil phase angle, (i.e., the location where a single 
hypothetical gradient coil would be located to produce an equivalent gradient pulse). 

Depending upon the geometric design of the gradient coils 76, 94, the duration, and the 
strength of the current pulse applied to the coil 76, 94, the spins in a sensitive region (e.g., one 
zone) will phase encode in one of the following manners: radial phase encoding, azimuthal phase 
encoding, axial phase encoding, or incomplete phase encoding. Firing a current of sufficient 
magnitude through the gradient coils causes the additional phase shift of the spins subjected to 
the magnetic field gradient of the coil to be distributed over a range from -180° to 180° and 
possibly vary by several multiples of 360° over the sensitive region. For radial, axial, and 
azimuthal phase encoding, a varying additional phase angle is generated so that the response of 
spins in a sensitive region (e.g., one quadrant) averages to zero. 

Azimuthal Phase Encoding (Dipole) 

Figure 18 illustrates a tool 90 with a dipole gradient coil 76 having a magnetic field 
pattern given by: 

B Gl i r > <P) = A cos <pr + A 9 sin <pq> 

where A r is the amplitude of the field along the radial (r) direction, Aq> is the amplitude along the 
azimuthal (9) direction, and (p is the azimuthal angle. The magnetic field gradient is oriented 
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circumferentially into the investigation region relative to the wellbore. The radial component of 
the field is responsible for phase-encoding the nuclear spins and is graphically illustrated in 
Figure 19, where the amplitude of the radial component is shown as a function of the azimuthal 
angle 9. 

The alteration of the nuclei spin rotation by the phase encoding gradient pulses is not 
equal across all of the azimuthal angles. The gradient magnetic field comprises a maximum 
amplitude at (p=0; reduces to approximately half amplitude at +60° and -60°; is approximately 
zero at +90° and -90°; is approximately negative one half amplitude at +120° and -120°; and is a 

o o 

negative maximum amplitude at +180 and -180 . 

For example, an NMR image utilizing five phase encoding gradient pulses can result in 
five bins of data, relating to various combinations of the eight zones of interest within the 
formation. The angular coverage of each zone can be different, since the gradient magnetic field 
strength varies as the cosine of the azimuthal angle. A graphical illustration of the zones with 
this particular phase encoding is shown in Figure 20. Two of the data bins will provide the spin 
densities directly for the top segment 120 and bottom segment 122, while the other bins will 
contain combinations of the spin densities of the other segments. 

The azimuthal resolution that can be obtained depends on the number of independent 
phase encoding gradient pulses utilized. The azimuthal resolution obtained from five phase 
encoding gradients provides eight zones as shown in Figure 20 and provides a top quadrant with 
an azimuthal angle of approximately 106° (+53° to -53°). Increased azimuthal resolution can be 
obtained with increased numbers of phase encoding gradients. 
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Radial Phase Encoding 

Figure 21 illustrates a gradient coil 76 with geometry useful for radial phase encoding. 
This is achieved by generating a strong gradient field with a single coil so that the additional 
phase due to the field of the gradient coil 76 varies within the thickness (i.e., shell) of the 
5 sensitive region by at least n 9 (wherein % radians is equivalent to 180° ). Figure 22 graphically 
depicts the magnetic field strength across the shell. Figure 23 illustrates a wellbore 12 and the 
partitioning of a section of formation 15 into a plurality of radially segmented sections 125a- 
125d. The magnetic gradient field is oriented radially into the investigation region relative to the 
O wellbore 12. In the portion of the shell close to the gradient coil, spins rotate faster than spins 
15 toward an outside portion of the shell. When the NMR signals are analyzed utilizing Fourier 
*! transforms, data can be grouped into a plurality of bins, the content of each bin representing the 

. 

measured signals from at least one of the radially segmented sections. The binned data can then 
|i be mapped to one or more radially segmented formation sections of the wellbore. 

tr: r 

15 Axial Phase Encoding 

Axial phase encoding can be used to obtain greater axial resolution. Figure 24 illustrates 
an example of gradient coil geometry useful for axial phase encoding. The phase angle varies 
along the longitudinal axis of the tool. Figure 25 depicts the magnetic field strength along the 
length of the tool. Figure 26 illustrates a wellbore 12 and the partitioning of a section of 
20 formation 15 into a plurality of axially segmented sections 126a-126d. The magnetic gradient 
field is oriented axially into the investigation region relative to the wellbore 12. When the NMR 
signals are analyzed utilizing Fourier transforms, data can be grouped into a plurality of bins, the 
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content of each bin representing the measured signals from at least one of the axially segmented 
sections. The binned data can then be mapped to one or more axially segmented formation 
sections of the wellbore. 

The present invention enables improved resolution imaging of a formation using a phase 
encoding means that produces altered NMR signals. These altered NMR signals can be analyzed 
utilizing Fourier transforms. Applying a pulsed magnetic field gradient within a particular 
direction, before the NMR signal acquisition, can alter the nuclei spins within the investigation 
region in that particular phase encoding direction. This gradient pulse is stepped through a series 
of differing amplitudes on subsequent excitations, so that the phase of the NMR signals along 
that direction are also incremented. Because of the position dependent nature of the NMR signal 
phase increments, the signals can be resolved by utilizing Fourier transforms along that direction. 
U.S. Pat. No. 6,326,784 assigned to Schlumberger Technology Corporation discloses one means 
of utilizing Fourier transforms of NMR signals that can be grouped within bins, the binned data 
representing groupings of NMR signals obtained from one or more azimuthal zones around the 
NMR tool. 

The NMR signal, /(G 9 ) for a particular applied gradient strength G is given by the 
following equation: 
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where /(q>) is the azimuthal formation spin density, to is the time duration of the phase encoding 
gradient, cp is the azimuthal angle, and y is the gyromagnetic ratio for the particular NMR 
sensitive nuclei. 

A term, k-space, is used in NMR imaging and it is defined as the following formula: 

The k-space defines the domain in which the raw magnetic resonance measurements are 
obtained. The amplitude of the phase encoding gradient is typically incremented with each 
excitation and correspondingly the phase of the NMR spins also increment. To satisfy the 
Nyquist requirement, the phase of spins is incremented by % radians in the imaged volume with 
each measurement. 

The NMR signal for a particular k<p value can be written as: 

The azimuthal formation spin density, which can also be referred to as the formation 
image /(cp), is the Fourier transform of /(k<p), which is written as: 

Kv)=jf(k 9 )e***dk 

The order in which the k-space is sampled is commonly either sequential or centric, 
although other sampling orders can also be used. In sequential ordering the phase encoding 
gradients are incremented linearly. In centric ordering, the phase encoding gradients are 
alternatively varied in the positive and negative direction. The signal to noise ratio of the NMR 
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image is largely determined by the low frequency components of the NMR signal (which are 
located close to a central position within the k-space). The high frequency components 
determine the sharpness of the NMR image (located at the peripheral region within the k-space). 

In one embodiment of the present invention, the method comprises applying a series of 
magnetic field gradients to a spin-echo series such that consecutive echoes are altered by a 
magnetic field gradient of differing amplitudes. The magnetic field gradient will effect the nuclei 
within the formation differently, depending on the nuclei position relative to the magnetic field 
gradient. The series of altered echoes can be detected and analyzed with respect to the differing 
magnetic field gradient amplitudes and segmented in relation to the azimuthal, radial or axial 
position relative to the wellbore. Utilizing this method, formation NMR responses can be 
determined for various segments within the formation. 

Embodiments of the invention can be utilized during the drilling of the wellbore for early 
formation evaluations, while other embodiments can be used in existing wellbores and deployed 
by wireline means. The invention can be used in conjunction with other instruments and tools, 
1 5 providing flexibility that can be used at the discretion of the operator. 

In further embodiments of the invention, multiple series of altered echoes, either from 
additional runs of a single gradient coil or the use of multiple gradient coils, can be generated, 
collected and analyzed. By analyzing multiple series of altered echoes, an image of the 
formation having a particular resolution (e.g., degree of angular segment) can be developed. 
20 When coupled with a MWD apparatus, the various segments of the formation can be translated 
into surface directional coordinates for use in formation evaluation and reservoir management. 
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Thus, the present invention enables the generation and imaging of NMR formation 
responses. Through use of the present invention, the resolution of the NMR imaging can be 
determined and controlled, thereby providing improved knowledge of the formation properties 
and assisting with an overall formation evaluation and reservoir management effort. The ability 
5 of the present invention to obtain formation properties without the use of radioactive sources 
reduces potential environmental concerns that accompany various other methods of formation 
evaluation. 



The foregoing descriptions of alternate embodiments of the present invention have been 
B presented for purposes of illustration and description. It is not intended to be exhaustive or limit 
f& the invention to the precise form disclosed. With the acquisition of formation imaging having 
improved resolution, the present invention may perform be useful in determining various 

# 

s formation properties, e.g., porosity, bound fluid volume (BFV), and permeability measurements. 
S It is also possible to perform azimuthal magnetic resonance imaging, which is useful for 
j£ interpreting heterogenous formations and performing geologically based steering in deviated or 

fa | 

ft horizontal boreholes. Obviously, many modifications and variations will be apparent to those 
skilled in the art. For example, the functionality of the MWD tool 52 may also be performed in 
the LWD tool 50 or divided between the MWD 52 and LWD 50 tools. Also, with a wireline 
tool, the gradient coil(s) may be located on a pad connected to the tool. Those skilled in the art 
will appreciate that the method and gradient coil(s) of the subject invention can be useful for 

20 eliminating the magnetic resonance signal of the borehole fluids, obtaining axially resolved NMR 

measurements, or NMR measurements with improved vertical resolution. For example, the 

length of the recessed area 70 along the longitudinal axis of the borehole can define the axial 

extent of an investigation region. A gradient coil or a plurality of gradient coils can be oriented, 

i 
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at known positions along the longitudinal axis of the borehole, within the recessed area 70. A 
current pulse applied to the gradient coil(s) will phase encode the spins in an axial segment of the 
formation. The embodiments were chosen and described in order to best explain the principles 
of the invention and its practical application, thereby enabling others skilled in the art to 
understand the invention for various embodiments and with various modifications as are suited to 
the particular use contemplated. It is intended that the scope of the invention be defined by the 
accompanying claims and their equivalents. 
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